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Abstract 

We present several searches for charmless hadronic two-body and three-body decays of B mesons 
from electron-positron annihilation data collected by the BABAR detector near the T{4S) resonance. 
We report the preliminary branching fractions B{B^ tt~^tt~) = (4.1 it 1.0 it 0.7) x 10~^, B{B^ 
K+TT-) = (16.7 ± 1.6 ± 1.3) X 10-6, B{B° p^n^) = (49 ± 13^^) x 10-^ B{B+ ij'K+) = 
(62 ± 18 ± 8) X 10-6, and present upper limits for nine other decays. 
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1 Introduction 



The study of B meson decays into charmless hadronic final states plays an important role in the 
understanding of CP violation. In the Standard Model, all CP-violating phenomena are a conse- 
quence of a single complex phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix 
|jl|. The Belle and BABAR collaborations have presented results [Q, ^ on measurements of CP- 
violating asymmetries in B decays into final states containing charmonium, leading to constraints 
on the angle P of the CKM unitarity triangle. Measurements of the rates and CP asymmetries for 
B decays into the charmless final states tttt and Kir can be used to constrain the angles a and 7 of 
the unitarity triangle Q. 

2 Detector and data 

Here we present new measurements of the branching fractions for charmless hadronic decays of 
B mesons in the final states^ 7t~^it~ and K~^7r^, and an upper limit for B K~ , which are 

based on a data sample consisting of an integrated luminosity of 20.6 fb^^ taken near the T(4S') 
resonance ( "on-resonance" ), corresponding to (22.57 ±0.36) x 10^ BB pairs. A data sample of 2.61 
fb""*^ taken at a center-of-mass (CM) energy 40 MeV below the T(4S') resonance ("off-resonance") 
is used for continuum background studies. 

We also report measurements of branching fractions and upper limits for B decays into non- 
resonant three-body modes and modes containing K* , p, u, and r]' resonances, which are based on 
a smaller data sample (see Sec. 

The data were collected with the BABAR detector at the PEP-II e^e^ collider at the Stanford 
Linear Accelerator Center. The collider is operated with asymmetric beam energies, producing 
a boost (/37 = 0.56) of the T(45) along the collision axis (z). For the analyses described in 
this Letter, the most significant effect of the boost relative to symmetric collider experiments is 
to increase the momentum range of the two-body B decay products from a narrow distribution 
centered at approximately 2.6 GeV/c to a broad, approximately fiat distribution extending from 
1.7 to 4.2 GeV/c. 

At this conference, the BABAR detector has been described in detail by Jim Panetta [||. Here 
we only reemphasize that the identification of tracks as pions or kaons is based on the Cherenkov 
angle 9c measured by a unique, internally refiecting Cherenkov ring imaging detector (DIRC). The 
K-TT separation varies as a function of momentum and is better than 8 standard deviations (cr) at 
1.7 GeV/c and decreases to 2.5a at 4 GeV/c. 

3 Analysis of Tr+Tr", K+tt", K+K' 

The selection of hadronic events for this analysis is based on track multiplicity and event topology. 
To reduce background from non-hadronic events, the ratio of Fox- Wolfram moments Q H2/HQ is 
required to be less than 0.95 and the sphericity 0] of the event is required to be greater than 0.01. 

All tracks are required to have a polar angle within the tracking fiducial region 0.41 < 9 < 2.54 
rad and a 9c measurement from the DIRC. The latter requirement is satisfied by 91% of the tracks 
in the described fiducial region. We require a minimum number of Cherenkov photons associated 
with each 9c measurement in order to improve the resolution. The efficiency of this requirement is 

^Charge conjugate states are assumed throughout, except where exphcitly noted. 



97% per track. Tracks with a O^. within 3cr of the expected value for a proton are rejected. Electrons 
are rejected based on specific ionization [dE/dx) in the DCH system, shower shape in the EMC, 
and the ratio of shower energy to track momentum. 

The kinematic constraints provided by the T(4S') initial state and relatively precise knowledge 
of the beam energies are exploited to efficiently identify B candidates. We define a beam-energy 



substituted mass tties = y ~ Ps' where E'b = (s/2 + pj • pB)/Ei, and ^/s and Ei are the 
total energies of the e~^e~ system in the CM and lab frames, respectively, and Pi and pb are the 
momentum vectors in the lab frame of the e~^e~ system and the B candidate, respectively. The ttt-es 
resolution is dominated by the beam energy spread and is approximately 2.5 MeV/c^. Candidates 
are selected in the range 5.2 < rriEs < 5.3 GeV/c^. 

We define an additional kinematic parameter AE as the difference between the energy of the 
B candidate and half the energy of the e~^e~ system, computed in the CM system, where the pion 
mass is assumed for all charged B decay products. The AE distribution is peaked near zero for 
modes with no charged kaons and shifted on average —45 MeV (—91 MeV) for modes with one (two) 
kaons, where the exact separation depends on the laboratory kaon momentum. The resolution on 
AE is about 26 MeV. Candidates with |Aii^| < 0.15 GeV are accepted. 

Detailed Monte Carlo simulation, off-resonance data, and events in on-resonance mss and AE 
sideband regions are used to study backgrounds. The contribution due to other S-meson decays, 
both from 6 — > c and charmless decays, is found to be negligible. The largest source of background 
is from random combinations of tracks and neutrals produced in the e~^e~ — > qq continuum (where 
q = u, d, s or c). In the CM frame this background typically exhibits a two-jet structure that 
can produce two high momentum, nearly back-to-back particles. In contrast, the low momentum 
and pseudoscalar nature of B mesons in the decay T{4S) — > BB leads to a more spherically 
symmetric event. We exploit this topology difference by making use of two event-shape quantities. 
The variable we considered that has the greatest discriminating power is the angle 9s between 
the sphericity axes evaluated in the CM frame, of the B candidate and the remaining tracks and 
photons in the event. The distribution of the absolute value of cos 6s is strongly peaked near 1 for 
continuum events and is approximately uniform for BB events. We require | cos^sl < 0.9, which 
rejects 66% of the background that remains at this stage of the analysis. 

The second quantity used in the analysis is a linear combination of the nine scalar sums of 
the momenta of all tracks and photons (excluding the B candidate decay products) flowing into 
10° polar angle intervals coaxial around the thrust axis of the B candidate, in the CM frame 
(Fisher discriminant Q Monte Carlo samples are used to obtain the values of the coefficients, 
which are chosen to maximize the statistical separation between signal and background events. No 
restrictions are placed on J^. Instead, it is used as an input variable in a maximum likelihood fit, 
described below. 

Signal yields are determined from an unbinned maximum likelihood fit that uses the following 
quantities: tues^ AE, T , and Qc- The likelihood for a given candidate j is obtained by summing the 
product of event yield and probability Vk over all possible signal and background hypotheses k. 
The rifc are determined by maximizing the extended likelihood function £: 





(1) 



The probabilities Vk(xj\OLk) are evaluated as the product of probability density functions (PDFs) 
for each of the independent variables Xj , given the set of parameters . Monte Carlo simulated data 



is used to validate the assumption that the fit variables are uncorrelated. The exponential factor in 
L accounts for Poisson fluctuations in the total number of observed events A^. For the X^vr^ terms 
the yields are rewritten in terms of the sum + and the asymmetry A= {n^ — nf)/{nj + nf), 
where nj {nj) is the fitted number of events in the mode B ^ f (B ^ f). 

The parameters for both signal and background ttt-es, AE, and PDFs are determined from 
data, and are cross-checked with the parameters derived from Monte Carlo simulation. The 6c PDFs 
are derived from kaon and pion tracks in the momentum range of interest from approximately 42 000 
D*^ D^TT^ {D^ K~TT~^) decays. This control sample is used to parameterize the 9c resolution 
as a function of track polar angle. 

The results of the fit are summarized in Table [l|. For the decay B^ — > K^K~ we measure 
the branching fraction B = (0.85t|jJ^ ± 0.37) x 10"^ The 90% confidence level upper limit for 

this mode is computed as the value for which Jq'' Ciaa^drik/ Cmaxduk = 0.90, where £max 
is the likelihood as a function of n^, maximized with respect to the remaining fit parameters. The 
result is then increased by the total systematic error. The reconstruction efficiency is reduced by its 
systematic uncertainty in calculating the branching fraction upper limit. The statistical significance 
of a given channel is determined by fixing the yield to zero, repeating the fit, and recording the 
square root of the change in — 21n£. 

Table 1: Summary of results for reconstruction efficiencies (e), fitted signal yields (Ns), statistical 
significances, and measured branching fractions {B). The total number of events entering the ML 
fit is 16032. For the K^K~ mode the 90% confidence level (CL) upper limit for the branching 
fraction is quoted. Equal branching fractions for T(4S') B^B^ and B^B^ are assumed. 

Decay mode e (%) Ng Stat. Sig. {a) B{m-^) 

B^^TT+TT- 45 41 ±10 ±7 4.7 4.1 ±1.0 ±0.7 

B^^K+TT- 45 169 ±17 ±13 15.8 16.7 ±1.6 ±1.3 

B^ K+K- 43 8.2l^;^ ± 3.5 1.3 < 2.5 



Event-counting analyses, based on the same variable set Xj as used in the fits, serve as cross- 
checks for the ML fit results. The variable ranges are generally chosen to be tighter in order to 
optimize the signal-to-background ratio, or upper limit, for the expected branching fractions. We 
count events in a rectangular signal region in the Ai^-m-ES plane, and estimate the background 
from a sideband area. The branching fractions measured using this technique are in good agreement 
with those arising from the ML fit analysis. Figure |l] shows the distributions in rriEs for events 
passing the tighter selection criteria of the event-counting analyses. 

The following sources of systematic uncertainty have been considered: imperfect knowledge of 
the PDF shapes, which translates into systematic uncertainties in both the fit yields and asymme- 
tries; systematic uncertainties in the detection efficiencies, which affect only the branching ratio 
measurements; and possible charge biases in either track reconstruction or particle identification, 
which affect only the asymmetries. 

The PDF shapes contribute the largest source of systematic uncertainty. Systematics due to 
uncertainties in the PDF parameterizations are estimated either by varying the PDF parameters 
within lo" of their measured uncertainties or by substituting alternative PDFs from independent 
control samples, and recording the variations in the fit results. The largest systematic uncertain- 
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Figure 1: The rngg distributions for candidates passing the selection criteria of the event-counting 
analyses for the tt+tt" (top), K^7r~ (center), and K^K~ (bottom) mode. 



ties of this type vary across decay modes and are in the range l%-7% . The total systematic 
uncertainties in the signal yields due to PDF systematics are given in Table |^. 

The overall systematic uncertainties on the branching fractions as given in Table Q are computed 
by adding in quadrature the PDF systematics and the systematic uncertainties on the efficiencies. 

4 Non-resonant three-body modes and quasi-two-body modes 

In Table |^ we summarize some earlier results that were obtained with a smaller datasample of 
7.7 fb^^ on-resonance, corresponding to 8.8 million BB pairs, and 1.2 fb^^ off-resonance. For the 
details of these event-counting analyses we refer to Ref. 9. 

Table 2: Summary of branching fraction measurements. Inequality denotes 90% CL upper limit, 
including systematic uncertainties. 
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5 Conclusion 



In summary, we have measured branching fractions for the rare charmless decays i? tt+vt , 
K^Tr~ , — > /C^TT^, and B^ — > r]'K~^ , and set upper hmits on B^ — > t:~^tt~t:^ , B^ P^t^~^ , 
B+ K+TT~7r+, B+ p^K+, B° j]' , B+ ujK+/tt+, B+ ioK°, and B+ K*^tt+. Our 
results are in good agreement with earlier measurements [ |10| . 
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